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ABSTRACT 

Aims. We present the main steps that will be taken to extract Ha emission flux from Javalambre Photometric Local Universe Survey 
(J-PLUS) photometric data. 

Methods. For galaxies with z < 0.015, the Ha-r[Nii] emission is covered by the J-PLUS narrow-band filter F660. We explore 
three different methods to extract the Ha -l- [N ii] flux from J-PLUS photometric data: a combination of a broad-band and a narrow- 
band filter (r' and F660), two broad-band and a narrow-band filter (/, i' and F660), and an SED-fitting based method using eight 
photometric points. To test these methodologies, we simulated J-PLUS data from a sample of 7511 SDSS spectra with measured Ha 
flux. Based on the same sample, we derive two empirical relations to correct the derived Ha-t[N ii] flux from dust extinction and [N ii] 
contamination. 

Results. We find that the only unbiased method is the SED fitting based method. The combination of two filters underestimates the 
measurements of the Ha + [N ii] flux by 22%, while the three filters method are underestimated by 9%. We study the error budget of 
the SED-fitting based method and find that, in addition to the photometric error, our measurements have a systematic uncertainty of 
4.3%. Several sources contribute to this uncertainty: the ditferences between our measurement procedure and that used to derive the 
spectroscopic values, the use of simple stellar populations as templates, and the intrinsic errors of the spectra, which were not taken 
into account. Apart from that, the empirical corrections for dust extinction and [Nii] contamination add an extra uncertainty of 14%. 
Conclusions. Given the J-PLUS photometric system, the best methodology to extract Ha + [N ii] flux is the SED-fitting based 
method. Using this method, we are able to recover reliable Ha fluxes for thousands of nearby galaxies in a robust and homogeneous 
way. Moreover, each stage of the process (emission line flux, dust extinction correction, and [N ii] decontamination) can be decoupled 
and improved in the future. This method ensures reliable Ha measurements for many studies of galaxy evolution, from the local star 
formation rate density, to 2D studies in spatially well-resolved galaxies or the study of environmental effects, up to m,./ =21.8 (AB; 
3<r detection of Ha-t[N ii] emission). 

Key words. Galaxies:interactions - Galaxies: statistics 


1. Introduction 

One of the most important processes driving the evolution of 
galaxies is the rate at which their gas is transformed into stars, 
namely the star formation rate (SFR). This parameter accounts 
for the amount of gas that is transformed into stars per unit time 
in a galaxy or a region of a galaxy. When the SFR is computed 
within a given cosmological volume, a SFR density (psfr) is 
obtained. 

The most direct indicator of the star-forming process is the 
ultraviolet (UV) radiation produced in the photosphere of young 
stars with intermediate to high masses (M > 3Mq). These stars 
emit more energy at short wavelengths (A < 3000A) than at 
longer wavelengths. As these stars are massive, their lifetimes 
are shorter than 300 Myr. Hence, UV traces star formation 
episodes of these timescales. The UV photons are, however, 
severely affected by dust extinction, which makes it difficult to 
translate UV flux into an SFR. Additionally, Earth’s atmosphere 
blocks most of UV radiation, which makes ground-based 
surveys in this wavelength range impossible. 

The UV radiation field of the most massive O stars 
(M > 20Mq) ionize the gas surrounding newborn stars. 
Hydrogen recombination leads to emission lines in different 
wavelength ranges, creating the Balmer series, among others. 
Thus, another widely used SFR indicator is the Ha emission that 


is originated in the nebular cloud surrounding the star-forming 
region. This emission is prominent during the first ~ 10 Myr of 
the star formation process, as the stars that can produce ionising 
photons die in these timescales. Hence, Ha emission probes 
more recent bursts than the UV. Moreover, this emission is less 
affected by dust than the UV continuum. Several calibrations 
exist to relate Ha flux with the SFR (e.g. |Kennicutt||1998[ 
|Calzetti|2013] l. 


Apart from these, there are other indirect SFR indicators. 
For example, the infrared (IR) thermal emission of dust heated 
by the UV field (see, for instance, Calzetti et al. 20071 or 


Perez-Gonzalez et al. 2008|l, or the forbidden {Oi\\AXii 2 i,V 29 


doublet ( [Kewley et al. |2004[ |Sobral et al.||2012) l. The latter is 
useful for studies in the optical for objects at z > 0.4, when Ha 
emission is shifted to the IR. 

The SFR has been studied looking at Ha for more than 
20 years, both in the local Universe and at higher redshifts. 
For example, the study of Gallego et al. (1995 1 computed 
the SFR density in the local Universe with the UCM Survey 
dat a ([Zamorano et al. 1994[ Gallego et al.| 1995| l. The work 
by Ly et al. ( |2007 1 uses narrow-band photometry to derive 
the luminosity function (LF) and the SFR at 0.07 < z < 1.47 
using information of Ha and [On]. Finally, the High Redshift 
Emission Line Survey (HiZEFS, Geach et al.|2008 1 also studied 
Ha emission at z = 2.23,1.47,0.84 and 0.40, using narrow-band 


1 



































Vilella Rojo et al.: Robust estimation of Ha flux with J-PLUS data 


photometry to estimate the SFR evolution of the last ~ 11 Gyr 
of the Universe ( |Sobral et al.|2013l l. 

With the Javalambre Photometric Local Universe Survey (J- 
PLUS0 Cenarro et al. in prep) we aim to derive the SFR of the 
local Universe measuring the Ha emission of nearby galaxies 
{z < 0.015). The design and survey strategy of J - PLUS allows 
us to probe the faint end of the LF(Ha). It is expected to reach 
~ 10^* erg ■ s ', which is 2.5 orders of magnitude deeper than 
the UCM Survey and the KPNO International Spectroscopic 
Survey (KISS, Salzer et al. 2000| 2001| l. In addition, thanks 
to the large surveyed volume, the bright end of the luminosity 
function is also constrained. 


The goal of this paper is to find the best method to ex¬ 
tract Ha fluxes from the J-PLUS photometric data. We present 
the methodology to analyse J-PLUS spectral energy distribu¬ 
tions (SEDs), and the estimation of the errors that arise with 
this methodology. The paper is organized as follows: Section]^ 
summarizes the J-PLUS photometric system and the main char¬ 
acteristics of the J-PLUS survey. Sectionpresents the dif¬ 
ferent methodologies that can be used to recover emission line 
fluxes from photometric data with a combination of broad-band 
and narrow-band filters, and in concrete the emission of Ha -H 
[Nii]. In Section]^ we test these methodologies with simulated 
J-PLUS photo-SEDs, based on a sample of 7511 Sloan Digital 
Sky Survey (SDSS) spectra of star-forming galaxies with mea¬ 
sured Ha fluxes. In Section|^we study the errors involved in the 
methodology. Dust and [N ii] corrections are treated in Section]^ 
Einally, in Section]^ we give the most relevant results and con¬ 
clusions of the work. Throughout this work, magnitudes are 
given in the AB system (Oke & Gunn| 19831. 


2. The J-PLUS Survey 


The J-PLUS survey is an auxiliary survey of the Javalambre- 
PAU Astrophysical Survey (J-PAS; Benitez et al. 20141. J-PLUS 
will be carried out at the Observatorio Astroffsico de Javalambre 
(OAJ; Cenarro et al. |2014|l using the Javalambre Auxiliary 


Surve y Telescope (JAST/T80) and TSOCam (Marin-Franch et al. 


2012 I. This is a wide-held camera (2.1 deg^ per exposure) in¬ 


stalled at the Cassegrain focus. It is equipped with a 9.2k-by- 
9.2k CCD with a pixel scale of 0.55"pixel '. 


Expected to start in 2015, J-PLUS will cover an area of 
~ 8500 deg^ of the northern sky, and has been particularly de¬ 
signed to carry out the photometric calibration of J-PAS using a 
specihc set of 12 hlters: 5 broad-band, and 7 narrow-band hlters 
located at key stellar spectral features. Pig. shows the trans¬ 
mission curves of the hlters. This hlter set allows us to retrieve 
accurate SEDs for millions of calibration stars, which will be 
used to transport the photometric calibration from J-PLUS to J- 
PAS dGruel et al.|20T2l l. 

In addition to the calibration goals, J-PLUS has been de¬ 
signed to acquire the Ha hux of the galaxies in the nearby 
Universe (z < 0.015) using one of the narrow-band hlters (F660) 
up to a detection magnitude ~ 22.5 in that band (see Table[T]for 
details of the J-PLUS photometric system). Moreover, the large 
held of view and efficiency of JAST/T80 and T80Cam make J- 
PLUS a powerful survey for other galaxy evolution studies, par¬ 
ticularly those that require spatially well-resolved galaxies. 


' http://www.j-plus.es 



Wavelength A 

Fig. 1. J-PLUS photometric system. Solid lines are the theo¬ 
retical transmission curves of the broad-band J-PLUS hlters. 
Coloured areas denote the theoretical transmission curves of 
the narrow-band J-PLUS hlters. The dashed line represents the 
quantum efficiency of the TSOCam. 


3. Methodologies 

In this section, we present a collection of methods that can be 
used to obtain the hux of Ha -i- [N ii] 4/16548, 6584 with J-PLUS 
photometric data. The width (~ 150A) and central wavelength 
(6600 A ) of F660 hlter ensures that we enclose the three lines 
inside the hlter at z < 0.015. We present the basic assumptions 
of each method and the most relevant equations in the following 
sections. 


3.1. Two filters (2F) method 


The simplest method that can be used is a combination of two 
hlters: one to trace the continuum, and another to contrast the 
emission line. This can be achieved with either two adjacent 
or overlapping hlters. Given the J-PLUS hlter system, we test 
the case that involves a broad hlter to trace the continuum (r’), 
and a narrow hlter placed at the wavelength range of the line 
of interest (7^660). This methodology is widely used in many 
photometric studies (e.g. |Ly et al.pOOTl [Takahashi et al.||2007 


Villar et al.|2008 Koyama et al.|2014 |An et al.|2014[); se e also 


the studies of the Ha galaxy survey (James et al. 2004| l for a 
combination of non-overlapping narrow-band hlters. 


In this case, and taking into account that r’ contains the hux 
of the Ha + [N ii] lines, the hux of these emissions can be recov¬ 
ered following the standard recipe (see Pascual et al.|2007' for a 
detailed description): 


(UFeeo - U,.-) 

F'ho'+[NII] = Af660-— > (1) 

where Ff^^o and F^' are the observed average Huxes inside the 
F66Q and r' hlters, and is dehned for any passband x at any 
wavelength of interest /Is as 

f Px (4) MA 
= paaTaJI, ’ 

where is the transmission of the passband x as a function of 
wavelength. In our case, /Is = Auo at z = 0, i.e. /Is = 6562.8 A. 
For J-PLUS P'660 and r’ hlters, we found = 125.3 A and 
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Table 1. J-PLUS photometric system properties. 


Filter 

name 

Central 

wavelength (nm) 

FWHM 

(nm) 

mag (AB) 
S/N=3 

mag (AB) 
S/N=50 

Comments 

u'j 

348.5 

50.8 

23.00 

19.69 

In common with J-PAS 

F378 

378.5 

16.8 

21.58 

18.01 

[0 ii]; in common with J-PAS 

F395 

395.0 

10.0 

21.60 

18.01 

Ca H-rK 

F410 

410.0 

20.0 

21.62 

18.03 

Hd 

F430 

430.0 

20.0 

21.65 

18.01 

G-band 

g' 

480.3 

140.9 

23.23 

19.92 

SDSS 

F515 

515.0 

20.0 

21.67 

18.01 

Mg 

r' 

625.4 

138.8 

23.26 

19.99 

SDSS 

F660 

660.0 

13.8 

22.64 

19.11 

Ha-t-iN ii]; in common with J-PAS 

i' 

766.8 

153.5 

22.31 

18.88 

SDSS 

F861 

861.0 

40.0 

21.48 

18.01 

Ca Triplet 

-c 

911.4 

140.9 

21.51 

18.19 

SDSS 


A,./ = 1419 A. The strongest assumption of this approximation is 
a flat continuum inside the r’ Alter. 


3.2. Three filters ('3Fj method 


To solve the problem of the flat continuum assumption of the 
2F method, we can use two Alters to trace a linear continuum, 
and a narTow-band Alter to contrast the line. There are different 


configurations for this method: in the work by Kennicutt & Kent 
( 1983| l no overlapping between the three Alters occurs, while the 
case of the narrow-band Alter overlapping two broad Alters is 
studied in |Pascual et al.| ( [2007| ). 

Because of the Alter conflguration of J-PLUS, we test the 
case in which the continuum is infen'ed using r’ and while the 
emission is inside F66Q. As in the 2F method, the r' flux also 
contains the lines’ flux, and it is risen from the stellar continuum. 
Because of this, we use an analytic formula to remove the Ha -H 
[N It] contribution inside the r' Alter, i.e. 


Hq'+[NII] - 




( 3 ) 


where F are the observed average fluxes, and a and p are defined 
as 


^ fA^F,(A)dA ^ A,,PAA = A,) 

f (A) AdA ’ fPx (A) AdA ’ 

where in our particular case, = Aua = 6562.8 A. These equa¬ 
tions are detailed in the Appendix [A| 

3.3. SED-fitting method 


where F^ is the observed flux for each of the x Alter with its 
error 6x, is the flux inside the x Alter of the j template, and k 
scales the templates at the magnitude of the galaxy that we are 
fitting. We estimate this scaling parameter for each j template by 
minimizing Eq.|^ i.e. 


kj = 


MW] 




( 6 ) 


Template models are simple stellar populations (SSPs) taken 
from Bruzual & Chariot ( |2003 BC03 hereafter), with 40 ages 
from 1 Myr to 13.75 Gyr in logarithmic bins, a Salpeter (1955 1 
initial mass function, and six metallicities ranging from 0.2Zo to 
2.5Zq. We extinguish these models with a Calzetti et al. ( |2000 1 
law from E(B-V) = 0 to E{B-V) = 1 in steps of 0.05. In the end 
we have 4200 templates at rest-frame. These are convolved with 
the J-PLUS photometric system properties using PySynphot|^ 
The template fluxes are then 


fT,Px(A)AdA 

Tx = —r -• ( 7 ) 

JPx{A)AdA 

The fitting is carried out without taking the flux of the 
F66Q Alter into account, which contains the Ha + [N ii] emis¬ 
sion fluxes. To derive the emission flux, we approximate it as a 
single line described as a Dirac’s delta function. This is called 
“the infinite thin line approximation” (see, for instance |Pascual 
|et al.|2007) l. If we split the total F66Q flux into two components 
(continuum and emission) and introduce it into Eq. we get 


Ff660 - 


f (■LHff+[NlI] + -Lcont) FF660 (^) 'f^A 

f Pf660 (A) AdA 


( 8 ) 


This method benefits from all the J-PLUS Alters to infer the 
emission flux after fitting the stellar continuum of the galaxy. 
To do that, we compare simulated J-PLUS observations with a 
set of template models, and for each pair observation-template, 
we compute the value of the x^ function, defined as 



( 5 ) 


where F'Ha+[Nii] is the emission flux, and Fcont the continuum. If 
we assume that FHa+[Nii] = FHa+[Nii]d (A - AhA, and split Eq.|^ 
in two integrals, we get 

Ff660 ~ Fp(560, cont T . (9) 

^F660 

^ A Python programming language adaptation of the widely used 
Synphot, developed hy the Space Telescope Science Institute (STSCI) 
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Fig. 2. Example of a SED fitting plotted over the SDSS spec¬ 
tra. Red asterisks show the original input (note the risen r’ dot). 
Green dots denote the photometric points after r’ decontamina¬ 
tion. These points overlap the red ones except for the r', because 
the r’ correction does not take them into account. Magenta solid 
line indicates the Her continuum using the 2E method. Cyan solid 
line shows the Ha continuum using the 3E method. Dashed blue 
line denotes the best-fitting template if no r' decontamination is 
done. Black solid line is the best-fitting template after r’ decon¬ 
tamination. 


Because the F66Q filter overlaps with r’, the Ha -H [Nii] 
fluxes are also inside the r' filter, which leads to 


Fr' = F, 


1 


Hq-+[Nii] 


( 10 ) 


To remove the Ha -i-[N ii] contribution from the r' filter, and 


obtain a more reliable r’ continuum, using Eq. 10 we decontam¬ 
inate the r’ flux by subtracting the F'hq.+[nii] flux that is inferred 
with the 3E method. 

We present an example of an SED fitting in Eig.[^ and com¬ 
pare the continuum that results from the equations of the 2E and 
3E methodologies. 


3.4. Measurements and error estimation 

To estimate the final Ha -H [N ii] flux and its error for each 
methodology, we perform 500 Monte Carlo runs. To do that, 
we perturb each passband x flux within a Gaussian distribution 
with jj. - F^ and cr - 6^ and apply each methodology. We keep 
a record of the inferred Ha -H [N ii] flux and perturb the origi¬ 
nal data again. In the end, we have an array of 500 Ha + [N ii] 
flux measurements for each method, noted as Fhq.+[nii]- Our final 
measurement for the Ha + [N ii] flux is the median of this array, 

(Fhq.+[nii]) = median (Fhq.+[nII]), (H) 

while the photometric error 5phot associated with this measure¬ 
ment is the median absolute deviation (MAD) of this array. 


|et al.|[2OT3| ) in the 10th Data Release (DR 10). We excluded 


Barion Oscillation Spectroscopic Survey (BOSS, Dawson et al. 


20131 galaxies, as its targets are luminous red galaxies at z > 0.2. 


We selected all the objects that were classified as star forming 
by the Portsmouth group, according to a BPT diagram criteria 
( Baldwin et al.|19^ . Thus, we do not expect a significant AGN 
contamination in the sample. We applied a spectroscopic red- 
shift cut at z < 0.02. This redshift cut is chosen to probe a vol¬ 
ume big enough to have a large number of galaxies with similar 
properties as the expected J-PLUS sources. After applying these 
criteria, we are left with ~ 12000 spectra. 

Erom this sample, we retain only the spectra with an Ha 
equivalent width (EWh^) 12 A > EWh^. This cut was done 
assuming that J-PLUS cannot resolve, with a precision of 3cr, 
EW/r 660 ^ 12 A because of the errors in the determination of 
magnitude of the zero point. Erom Pascual et al. ( 2007| l we know 
that 


EW = Af66o(e-l)^^—(13) 
1 - Qe 

where e = Af66o/A,,, and 2.5 log Q = m^' -m/reeo- Assuming a sys¬ 
tematic error in the determination of the zero-point magnitude of 
dm 0.02 in r’ and F66Q, there is a limiting difference in mag¬ 
nitudes of dm ~ 0.03, which we cannot resolve with enough 
confidence, and which leads to a minimum EW of detection that 
is EW F 660 - 12 A. Erom the SDSS data we cannot know the ob¬ 
served EWFgeo, so we take EWh^ = 12 A as lower limit. Einally, 
a last cut in the median signal-to-noise ratio (S/N) of each spec¬ 
trum was done. Only spectra with average S/N > 5 are kept. 
This limit is chosen to guarantee that spectroscopic flux errors 
are not important. 

After applying these criteria, the sample contains 7511 spec¬ 
tra. We refer to this sample as SI. Eigure [^presents the distribu¬ 
tion of the SI as a function of the EW /?660 and m^/. As, for now, 
the redshift is not taken into account in our analysis, we shifted 
all the spectra of S1 to z = 0. After this, we convolved them with 
the J-PLUS filters that are in the wavelength range of the SDSS 
spectra to obtain F^. Because these spectra have a shorter wave¬ 
length range than J-PLUS filters, we lose the information of u'j, 
F378, and z' bands. The convolution retrieves the apparent mag¬ 
nitudes and mean fluxes of each passband. We stress that these 
apparent magnitudes are computed from the flux enclosed inside 
the fibre, and are not representative of the whole galaxy. With 
the apparent magnitudes, we compute the expected S/N using 
the J-PLUS exposure time calculatotj^ This tool provides the es¬ 
timated S/N given an apparent magnitude, an exposure time (that 
can be divided into several exposures) and different sky condi¬ 
tions. The assumed conditions to compute the S/N values were a 
grey night, a seeing of 0.9 arcsec, a photometric aperture of 1 .8 
arcsec, and an airmass of 1.2. With this, the flux error in each 
passband xisd^ - ^Vs/n,- 

We apply each of the methodologies explained in Sect.j^ to 
the J-PLUS photo-spectra of SI and study the performance of 
each methodology below. 


5phot - 1.48 X median (| FHa+[Nii] - (Fhq-+[nii]) I) ■ (12) 

4. Testing the methodologies 

4.1. Datasampie 

To test each methodology, we use a set of SDSS spectra with 
emission lines measured by the Portsmouth Group (|Thomas| 


4.2. J-PLUS vs. SDSS 

At this point we test the precision of each method described in 
Sect. I^to recover the flux of Ha -H [Nii]. To this aim, we com¬ 
pare our inferred fluxes with those provided by the Portsmouth 
Group. As their values are dust corrected, but the SDSS spectra 

3 www.cefca.es/jplusetc 
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Fig. 3. SDSS star-forming sources with z < 0.02 (grey areas); 
Top panel; red line shows the sample S1 distributed as a func¬ 
tion of mr'. Dashed area denotes the sample S2 (see Sect. 5.3. l| l 
distributed as a function of m,’. Bottom panel: red line shows 
sample SI distributed as a function of logEWna. Dashed area il¬ 
lustrates the sample S2 (see Sect. 5.3. 1[ ) distributed as a function 
of logEWHa. 


are not, we add dust to the SDSS measurements following the 
[Calzetti et al.| ( |20()0| l extinction law with Ry - 4.05, which the 
Portsmouth group applied, and the values of E(B - V) provided 
by them. 

We apply each methodology to the spectra of SI, and com¬ 
pare the recovered flux of Ha -i- [N ii] lines with the spectro¬ 
scopic measurements. 


P _ (Ehq.+[Nii]) 
^ “ rSDSS 

Ha+tNlI] 


where, in this case, is 

’ ’ Hq'+[NiiJ 


(14) 


Fig. 4. Normalized distribution of the ratios between the recov¬ 
ered and the spectroscopically measured Ha -i- [N ii] flux for S1 
spectra, with the 2F (empty red histogram) and 3F (dashed blue 
histogram) methods. Blue curve shows the Gaussian htting to 
the distribution of R in the 3F method. The best htting values are 
labelled in the panel. 


that is given by the broad hlter. We lose the information about 
the colour of the galaxy, i.e. the true shape of the continuum at 
the wavelength range of the emission line. The distribution of 
results is very asymmetric, and cannot be htted to a Gaussian 
distribution. The median of this distribution is yu = 0.78, while 
the dispersion is given by (T 2 f = 0.5 x (7*84 - Pie) = 0.20, where 
^84 and Pi6 are the 84* and 16* percentiles, respectively. 

To cope with this bias, |Sobral et al. ( |2009| 2012} , and subse¬ 
quent works from the HiZELS survey, and also the work by 
|et al.| ( |201 l| l, introduce a correction based on the broad-band to 
narrow-band colour. This correction is applied to compensate the 
flat continuum assumption and the difference between the central 
wavelength of the broad-band Alter and the narrow-band cen¬ 
tral wavelength. This correction would be a step between the 2E 
method and the 3E method. However, we prefer the 3E method, 
as its correction is analytic. 

With the 3E method, the distribution of results becomes al¬ 
most Gaussian, but a bias of ~ 9% still persists, with yu = 0.91 
and a dispersion cr 3 F = 0.06 from the Gaussian fit. In this case, 
we are more sensitive to the true shape of the continuum in the 
wavelength range of the line, but we still have to assume it is 
linear, which is a poor approximation because of the Ha absorp¬ 
tion. 


4.2.2. SED fitting 


■^Hff+fNii] “ + E[nII],46548 + E[n ii],46583] R 10 D 

(15) 

because we are adding dust attenuation to its measurements. The 
resulting distribution of ratios R is fitted to a Gaussian in all the 
cases. 

4.2.1. 2F and 3F methods 

We show the resulting distribution of R when we analyse the 
spectra applying the 2E and 3E methodologies in Eig. The 2E 
method is biased, mostly because we assume a flat continuum 


We plot the distribution of results after applying the SED fitting 
methods in Eig. We compare the results without applying the 
r’ decontamination. The SED fitting routine with r' decontami¬ 
nation performs better than the SED independent methodologies 
(i.e. the 2E and 3E methods), being unbiased (yUsEO = 1.00). 
With this technique, we do not approximate the continuum to 
any function; we use the continuum inferred from BC03 tem¬ 
plates, which also contain an estimation for the absorption of 
Ha . Another interesting result is the impact that has the Ha + 
[Nil] emission inside the r’ Alter. Not taking this into account 
causes a bias in our results by ~ 8%. 
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Fig. 5. Normalized distribution of the ratios between the recov¬ 
ered and the spectroscopically measured Ha -H [Nii] flux for 
SI spectra, with the SED fitting routines. Empty red histogram 
shows the results after SED fitting method without r’ decontam¬ 
ination. Dashed blue histogram denotes the results after SED 
fitting method with r’ decontamination. Solid curves indicate 
Gaussian fits to the distribution of R. The best fitting values are 
labelled in the panel. 


4.3. Testing the methodologies: conclusions 


At this point we conclude that the methodologies that only 
use two or three filters are inconvenient for J-PLUS, given our 
filter configuration, as they produce biased results in average. 
In the case of the 2E methodology with no colour correction, 
this bias is ~ 22%, while the 3E method underestimates the 
Ha -I- [N It] flux by 9%. As |Sobral et al. ( 2012| l points out, this 
effect when using only two filters is, in part, caused because the 
central wavelength of the narrow filter does not coincide with 
the broad-band central wavelength. To cope with this, we need 
more information, and so we include the i' filter. However, this 
does not solve totally the problem. 

With the SED fitting procedure we avoid this bias, though 
attention must be paid to the contribution of Ha -H [N ii] inside 
the r ’ filter. Not taking this effect into account biases our results 
by 8%. Einally, the SED fitting procedure with the r’ decontam¬ 
ination is unbiased and has a dispersion of ctsed = 0.06. This 
dispersion is a combination of the errors associated with the 
photometry, dphou and other factors that are discussed in detail 
in Sec. 15.41 


Fig. 6. Green dots denote the medians of the Gaussian fitting 
to the galaxies inside the magnitude bin which median is m^s 
Green bars indicate cr of this Gaussian fit; The shaded area is the 
4.3% uncertainty defined by 5syst- 


Figure shows the g. of these distributions as a function of 
the median m^- magnitude of each bin. Error bars are the stan¬ 
dard deviation cr of each fit. The results are well recovered in 
all the magnitude ranges. Error bars increase from cr ~ 5% at 
nir' - 15 to cr ~ 8% at mr' = 19.5. We interpret this dispersion 
as the combination of two effects. We see a more detailed study 


of these dispersions in Sect. 5.3 and Sect. 5.4 


5.2. Dependence on EW 

We repeat the same analysis, but binning our S1 as a function of 
the observed EWfceo- This EW is not the one that we used in the 
data selection criteria, but that resulting from the recovered flux 
ofHa-H [Nil], 


EWe 660 - 


{Fhq.+[Nii]) 

FF 660 - {F Hq-+[Nii]) 


(16) 


Results are shown in Fig.|^ Each bin contains 375 spectra. In 
this case, the error weighted median is still unbiased, though the 
results show a trend in the recovered flux that creates an excess 
in the region of small EW, and an underestimation in the larger 
EWs. This latter effect might be due to a bad determination of 
the continuum in regions dominated by ionized gas, where SSP 
models are not valid. Nevertheless, the error is constrained to 
under 4.3%. 


5. SED fitting routine: performance and error 
budget 

In the previous section, we conclude that the SED fitting rou¬ 
tine is the most reliable methodology for our purposes given the 
J-PLUS filter configuration. In this section, we perform some ad¬ 
ditional tests on this methodology and explore the error budget 
in the measurements. 


5.3. Simulating observations at highernir' 

We have shown that observing at different does not introduce 
any biases. Our spectroscopic sample however cannot reach the 
expected J-PLUS mr' limiting magnitude (the selection criteria 
naturally cut the sample before mr' = 20.5, as seen in Fig.[^. To 
study the performance of the SED fitting method at magnitudes 
higher than mj. = 20.5, we have to simulate observations. 


5.1. Dependence on m,' 

To see if any bias appears at faint magnitudes, we study the dis¬ 
tribution of R at several mr' bins. To do that, we selected the 
galaxies in SI within a mr' range and fit a Gaussian to their 
R distribution. Bins are defined to contain the same number of 
galaxies, which in this case is 626. 


5.3.1. Sample S2 

We selected a subsample of galaxies with mr' < 17.5 and me¬ 
dian S/N > 20 from the parent sample , which were artificially 
scaled at any magnitude. These thresholds were chosen to have 
a subsample of galaxies with good quality data and a reason¬ 
able number of spectra to do statistics. In particular, we selected 
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Fig. 7. Green dots represent the medians of the Gaussian htting 
to the galaxies inside the EW bin which median is EW; Green 
bars denote the cr of the Gaussian htting; Shaded area shows the 



15 16 17 18 19 20 21 22 

m,.' 


Fig. 8. Results after simulating observations of S2 spectra. Green 
dots indicate the medians of the Gaussian htting to the galaxies; 
Green bars denote the cr of this Gaussian htting; Shaded area 
illustrates 4.3% uncertainty dehned by (Jsyst- 


5.4. Estimating the errors 

Here we carry out an analysis of the errors associated with 
the measurement method. For convenience, we refer to the dis¬ 
persion of the results, given by the standard deviation of the 
Gaussian hts, with letter cr; we reserve letter h for the errors 
and uncertainties in the measurements. In Sect. |3.4| we explained 
how we estimate the photometric errors hphot of {Fhq.+[nii])- 
However, in Fig. we see that the dispersion of the results in 
the brightest magnitudes is almost constant. This means that we 
cannot explain the dispersion of the results only with dphot, and 
that we need to add a new uncertainty to our measurements, i.e. 

O-SED = + ^syst • (17) 


To compute the value of dsyst, we want to minimize 5phot- We 
simulate observations with S/N > 10^ for S2. In this case, the 
dispersion in the results is only due to 6syst- Doing this, we hnd 
that (5syst = 0.05. The shaded area in Figs. andshows this 
systematic uncertainty. We see that this value well constrains the 
error bars until these begin to increase because of photometric 
errors. 

To test the validity of this result, we compare this new er¬ 
ror with the purely photometric error. To do that, following the 
Monte Carlo approach described in Sect. 3.4 we create 1134 
normal distributions (one per spectrum in S2), centred in yu, = 1 
and with cr,- = dphot, i- We see in Fig.lOlthat, when we add all these 
distributions, we recover a new distribution whose dispersion is 
only due to photometric errors (blue solid line). However, if we 

repeat the same exercise with cr,- = -^^phot, i ^sysc 

the resulting distribution has a dispersion that resembles that of 

our results (red dashed line). 


Finally, we apply this systematic error to the 7511 spectra of 
SI. We show in Fig.[^the same distribution as in Fig.|5 but we 
overplot the distribution when only taking dphot (blue so id line) 
into account and when we add the 4.3% uncertainty. This new 
distribution (red solid line) traces the dispersion of our results 
well, meaning that our error budget is reliable. 


high S/N galaxies because we want J-PLUS photometric errors 
to dominate over the specttoscopic errors. We end up with a sub¬ 
sample of 1334 galaxies, which we call Sample S2 (Fig.[^. 

5.3.2. Simulation routine 

The spectra of S2 are scaled to any m,./ of interest. After we 
change the magnitude of each hlter with the same difference 
of magnitudes that we apply to match the inttinsic m^- with 
the desired mr', we compute the expected S/N for each filter, 
and we perturb each of the x fluxes within its error bar from 
a random normal distribution with ji - and cr = ^‘/s/n,. 
These perturbed fluxes are now considered the observations, 
and from these fluxes we estimate (fHa +iNni ) and its uncertainty 
following the process described in Sect. |3.4| 

Figure 1^ shows that we recover the fluxes without biases up 
to magnitude m,./ ~ 21.8. The errors increase at fainter magni¬ 
tudes, while remains constrained at a ~ 4.3% at brighter magni¬ 
tudes. We discuss this in the next section. 


5.5. SED fitting routine: performance and error budget 
conctusions 

In this section, we have studied the sources of error that affect 
our measurements of Her -H [N ii] flux when using the SED fit¬ 
ting method. We have seen that our measurements are not biased 
at any magnitude of the J-PLUS detection magnitude range. To 
accomplish this, we studied S1 spectra at their magnitudes and 
simulated observations of S2 galaxies at magnitudes in which we 
have no data. With this test, we And that there is an uncertainty 
in all the magnitudes that is independent of the photometric er¬ 
rors. 

We studied this uncertainty and treated it as a source of error. 
We see that adding a systematic uncertainty of 4.3% to the (5phot 
of each measurement, allows us to recover a distribution that re¬ 
produces the observed dispersion in the results at any magnitude. 
This uncertainty can be due to a combination of several sources, 
such as not taking the intrinsic errors of the SDSS spectra into 
account, fitting SSPs to regions that may not be well represented 
by SSPs, or differences in the measurement procedures between 
our results and the Portsmouth Group. 
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R 



R 



Fig. 9. Top panel: black histogram denotes the distribution of R 
for the 1134 galaxies in S2, scaled to m^/ = 17.5. Blue curve 
shows the sum of 1134 normal distributions centred at /t = 1 and 
with cr = dphot- Red dashed curve is the sum of 1134 normal dis¬ 
tributions centred at /t = 1 and with cr - Bottom 

panel: same as above, except for spectra scaled to m^/ = 21.25. 


6. Dust correction and [N ii] removal 

6.1. Dust correction 


Our aim is to recover the Ha emission from galaxies; however, 
galaxies contain dust that is mixed with the stellar populations 
and the gas. Dust is present in molecular clouds before they col¬ 
lapse to form stars, and it is mixed with the hot gas after the 
first stars of a star-forming region are born. The presence of dust 
has two important observable consequences: it attenuates the to¬ 
tal amount of light that we receive, and tends to redden the true 
colour of light-emitting region. Both effects have an impact on 
the absolute and relative fluxes and magnitudes that we measure. 
Attenuation causes all the magnitudes in the optical and UV to 
increase, while reddening causes that this increase is higher as 
we move to the blue parts of the spectrum. 

The proportion between the attenuation in the Johnson V 
band (Ay) and the difference between the observed and the in¬ 
trinsic (i.e. dust-free) magnitudes of the Johnson B and V bands 
(i.e. the colour excess E(B - V)) has been called the extinction 
law (see Cardelli et al.|1989] l, although definitions in other bands 
exist (Fitzpatrick|1999|l, i.e. 




Ay 

E(B - V) ■ 


(18) 


Fig. 10. Distribution of R for all the 7511 galaxies in SI. 
Blue curve denotes the sum of 7511 normal distributions cen¬ 
tred at ju = 1 and cr = dphot- Red dashed curve indicates the 
sum of 7511 normal distributions centred at yu = 1 and with 

V^phot + 


From [Calzetti et al.| ( |2000[ ), we know that the relation be¬ 
tween the intrinsic flux and the observed one is 

Fi(A) = Fo(/l)10°'‘'^-‘ = (19) 

where k'(A) is a polynomial that depends on A and Ry. 

To compare the recovered flux with the one provided by 
the Portsmouth Group, we added dust to their measurements 
with the values of E(B - V) that they provided (see Eq. 15 ). 


However, with the filter configuration of J-PLUS it is more 
difficult to estimate the dust contribution to our fluxes. 

From now on we will be using the SED fitting procedure 
to extract the flux of Ha -H [Nii], as the other methodologies 
presented biases. We stress that the corrections presented in 
this section can be modified in the future because they are 
independent of the method that we used to isolate the emission 
lines. 


A common assumption that is applied to photometric data 
is an attenuation of Ah^ = 1 mag for the Ha emission (see 
Kennicutt|[T992l IGeach et al.|[2008l IVillar et al.|[2008l |Sobral| 


etal.|2012 


_ . The dust-corrected values assuming this attenuation 

are RHa+[Nii] = Rfeeo x = 2.5Rf660- However, we ap¬ 

plied this correction and found that it tends to overestimate the 
results, being the median of R - 2.3, far from the expected 

R=l- _ _ 

Other studies, such as |Sobral et al.| ( [20T4) l, use more sophis¬ 
ticated techniques to correct dust extinction. In their work, the 
correction by Garn & Best (2010 1 for star-forming galaxies is ap¬ 
plied. This correction relates dust extinction with the stellar mass 
of galaxies. However, we do not aim to derive stellar masses for 
our test sample, as at the redshift range of our interest the SDSS 
spectra generally do not cover the entire galaxies, but a small 
portion of them is covered by the SDSS fibre. 

With S1 we find that there is a trend between the observed 
g' - i' colour and the spectroscopically measured E(B- V). This 
is represented in Eig. [m To fit a power-law function to these 
points, we bin this sample in g' - i' logarithmically spaced bins 
and compute the median of each bin. We see in Eig.[^that, for 
g' - i' < 0.5, medians of E(B-V) are 0. This is because there is 
a subsample of data, which has E{B - V) = 0 and no associated 
error. We take them into account as zero for the fit and we obtain 
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Fig. 11. Distribution of spectroscopically derived E(B - V) and 
observed g' - i'. Blue dots denote the medians of E{B - V) val¬ 
ues enclosed in the logarithmically spaced bins. Red dashed line 
shows the ht to the blue dots. 



Fig. 12. Relation between the spectroscopically measured and 
dust connected Ha flux and the total Ha -H [N ii] flux. We see that 
the distribution is bimodal. The two trends can be differentiated 
if we split the sample by its observed g' - i' colour. 


E(B - y) = 0.206 {g' - i') - 0.0457. 


( 20 ) 


We avoid negative colour excess values by making 
E(B - y) = 0 if the inferred value is negative. The derived 
relation should only be used in certain circumstances. The g' - i' 
colour that we are using here is obtained after convolving the 
SDSS spectra with the J -PLUS photometric system. This means 
that this colour is a local property of the region inside the SDSS 
hbre. With J-PLUS we will have spatially resolved galaxies, 
where we will be able to differentiate and isolate star-forming 
regions. These are the regions where this relation is reliable. For 
galaxies that are not spatially resolved, the g' - i' colour is an 
overall colour resulting of the underlying stellar populations and 
the gas, if any. In these cases, the validity of this correction is not 
ensured, and other SED-htting codes which study galaxy prop¬ 
erties and stellar populations in more detail, may be used to ex¬ 
plore dust extinction (see, for instance, MUFFIT, |Dfaz-Garcra 
let al.|2015» . 


6.2. /N iiy correction 

The F66Q hlter contains the flux of Ha and [N ii] doublet, and it 
is not possible to deblend these three lines to isolate the Ha emis¬ 
sion flux. To cope with this problem, empirical relations must be 
applied. With data from SI, we hnd that there is a relation be¬ 
tween the flux of Ha -H [N ii] and the Ha flux alone. Figure [T3] 
shows this relation. We hnd that there is a slight bimodality in 
the distribution of huxes, which is blurred in the low-emission 
regime. This bi-modality can be disentangled with the help of the 
colour g' - i', and we ht the following equation to each branch: 


log(Ha) = 


0.989 log(FHff+[Nii],D.c.) - 0.193, 
0.954 log(FHa+[Nii],D.c.) - 0.753, 


if g' - i' < 0.5, 
if g' - i' > 0.5, 
( 21 ) 


where F'Hff+[Nii], d.c. refers to the F660 Hux after dust correction. 


6.3. Ha only measurements 

We apply both corrections to our measurements and compare 
the recovered results with the spectroscopic values of Ha -i-[N ii] 
without dust and Ha only. As we can see in Fig. [T^ both cor¬ 
rections help us to recover the Ha Hux without adding any extra 



1.0 

R 


Fig. 13. The empty red histogram is distribution of the recovered 
Ha -H [N ii] after correcting for dust; dashed blue histogram de¬ 
notes the distribution of the recovered Ha hux after correcting 
for dust and [N ii]; Solid lines indicate Gaussian hts to the data. 
Best-htting values are labelled in the panel. 


bias to the whole set of measurements. To these results, we ht a 
Gaussian distribution. 

It is important to stress that [N ii] correction is empirical, and 
can only be applied after correcting for dust the observed F660 
hux. This relation should hold regardless of the dust correction 
that is applied, as it has been calibrated with dust-free data. 

6.4. Error budget 

After correcting the Hux of Ha + [N ii] from dust, the resulting 
distribution has a larger dispersion, although it does not become 
biased. This increase in the dispersion must be taken as another 
source of uncertainty, namely d^on, i-C- 

= -y/^phot + ^?yst + ^corr ■ (22) 

To derive the value of 5corr, we compare the dispersion in 
distribution of raw Ha + [N ii] (crgEo = 0.06, Fig. |^, only with 
the dispersion of the distribution of Ha (crHa = 0.15, Fig. [T3] l. 
We hnd that 5corr = 0.14. 


9 






















Vilella Rojo et al.: Robust estimation of Ha flux with J-PLUS data 



R 

Fig. 14. The Histogram is the distribution of the recovered Ha 
flux after correcting for dust and [Nii]; The solid curve shows 
the distribution only when taking (5phot into account; the dashed 
curve shows the distribution when adding dsyst and (Jcoit to dphot- 


galaxies by ~ 20%. The asymmetry of the distribution makes 
it difficult to treat it statistically and to have the errors under 
control. To use this methodology, colour corrections should 
be applied, as shown by |Sobral et al.| ( |2009| ). 

2. Using a combination of two broad-band filters and a narrow- 
band gives better results. However the method is still biased 
in average by a 9%. 

3. Fitting the whole SED to a collection of SSP models after 
subtracting the Ha -i- [Nii] flux from r’, retrieves unbiased 
Ha -I- [N It] fluxes. We stress the importance of taking into 
account the contribution of the emission lines inside the r’ 
filter. Not correcting the flux of this filter introduces a bias of 
8 %. 

We conclude that the SED fitting method is the best one 
given the J-PLUS capabilities. We encourage other photometric 
surveys targeting emission lines with narrow-band filters to 
explore the SED fitting methodology, instead of restricting them 
to the use of two or three filters. 


Combining both uncertainties dsyst (Sect. 5.4 1 and da 
obtain that the final error in Ha is 


we 


dviry — 


V^phot + 0.15U 


(23) 


We now add this error to each measurement of S1 spectra 
and repeat the same analysis as in Sect. 5.4 The resulting dis¬ 
tribution, with the fitting properties and the errors, is shown in 
Eig. 14 We see that a 15% uncertainty creates a distribution that 


resembles the dispersion of R after correcting for dust and [N ii], 
as desired. 


6.5. Dust correction and [N w] removal: conclusions 


To correct the observed Ha + [N ii] emission flux from dust 
and remove the [N ii] contribution, we derived empirical correc¬ 
tions from the SDSS data. After that, the recovered Ha flux is 
still unbiased, but suffers a larger dispersion. 

Einally, we demonstrate that the error of our measurements 
of Ha flux, diia, has three contributions. The first is given by the 
photometric errors of our data, dphot- The second contribution 
has several sources, which include not considering the intrinsic 
errors of SDSS spectra, the use of SSPs to fit the stellar contin¬ 
uum of regions that may not resemble SSPs, and differences be¬ 
tween our measurement procedure and the Portsmouth group’s 
method. This results in a systematic uncertainty, (5syst. The last 
source of uncertainty is related with the corrections of dust and 
[N It], namely 5corr- In the end, we express the error of our Ha 
measurements as 


In this section, we have studied how to correct for dust extinc¬ 
tion and for the contribution of the [N ii] doublet to the total ob¬ 
served flux inside F66Q. We used two empirical relations derived 
from SDSS data. Taking the properties of galaxies used to derive 
these expressions into account, they should be representative of 
the properties of star-forming regions of galaxies in the local 
Universe. 

When applied, both corrections retrieve unbiased results. 
This means that we can decouple them from our first goal: ob¬ 
taining reliable measurements of Ha + [N ii]. In this sense, both 
corrections can be modified in the future if better corrections are 
found. 

When we add both corrections to the raw measurement of 
Ha -H [N It], the distribution of results is unbiased, but the disper¬ 
sion increases. We obtain that the uncertainty introduced when 
we apply our corrections is ^con ~ 14%. 


7. Summary and conclusions 

We have presented the capabilities of the J-PLUS survey to infer 
Ha emission from photometric data. We first presented different 
methodologies and equations that can be applied to extract emis¬ 
sion line fluxes from narrow-band and broad-band photometry. 
After that, we tested each methodology simulating observations 
of SDSS spectra as seen by J-PLUS. We find that: 

1. Using a broad- and a narrow-band filter without taking care 
of the colour of the galaxy retrieves severely biased results, 
tending to underestimate the Ha -H [N ii] flux of star-forming 


^Hq- - 


^phot ^syst + <5corr> 


(24) 


where 5syst = 0.05 (Sect. 5.4 1 and dcon- = 0.14 (Sect. 6.4 1 . 


This means that our Ha measurements have a 15% of uncer¬ 
tainty, regardless of the quality of the data. 

We stress here that this 15% is an upper limit to the uncer¬ 
tainty related to the methodology and the corrections. Because 
each stage of the process (emission detection, dust correction, 
and [N n] correction) is unbiased, we can decouple them. This 
would allow us to improve these corrections and reduce the sys¬ 
tematic error. Eor instance, J-PLUS has a narrow filter in the 
[On] wavelength range. This forbidden emission is also a tracer 
of the SER, and could help us constrain the dust contribution 
better than the colour that we used. 
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Appendix A: 3F method: Equations 

Here we develo p th e equations for the 3F method that was de¬ 
scribed in Sect. 3.2[ The average flux integrated any filter x 
can be obtained, if the passband P properties are known, with the 
following expression: 


_ ^F,PAA)MA 

f = ± _ 

f PAX) AAA 


(A.l) 


where P^ is the transmission of the passband x, as a function of 
wavelength. 

Inside the r’ filter, there are two main contributions: the 
flux of the continuum Fr'^com, and the flux of Ha -i- [Nii], 
T^r',Ha+[Nii]- Following the same reasoning, in the FbbO filter 
there is Tfeeo , cont and F'f660,Ha+[Nii]- 

We approximate the continuum to a linear function 

i"contU) = Md + (V, (A.2) 

and for the emission, we assume that the flux of the three lines 
can be understood as one single, infinitely thin line centred at the 
Ha wavelength, which is described as a Dirac’s delta function 
(the so-called infinite thin line approximation): 

■Fline = ■FHff+[Nll]<J(d - /Ihq.). (A.3) 

With this, the observed flux inside r’: 


- f contP.' (d) AAA f PlinePF (d) ddd 

F I = - -h - - (A 4) 

f Pr' (d) AAA f Pr> (A) AAA 


Plugging Eqs. |A.2| and A. 3 into Eq. |A.4| and rewriting the 
integrals, we get 


_ f (MA + N) Pr' (d) AAA 

Fr> - - - -—- +Pr’F\me - +/3r'Flme + A , 

f Pr' (d) AAA 

(A.5) 

where a and j3 can be defined at any passband x, at any wave¬ 
length of interest dj, as 

^ /d^p,(d)dd ^ a.paa^a,) 

^ PAX) AAA JpAA)AAA' ’’ 

In our case, = dn^. Eollowing the same steps for the F66Q 
and the i’ filters, we obtain 


Ff 660 - Map(,(,0 + /?F660Pline + A , 


(A.7) 


F i' - Moi' + N . 


(A.8) 


In Eq. A.8 the line contribution does not appear because the 


i’ filter does not cover it. Combining Eqs. |A.7| and |A.8| we obtain 
Ff660 - F -,' - ^F660Pline 


M = 


<TF660 - O',/ 


Plugging this into Eq. A.8 
N = Fi' - a, 


FF 660 - Fi' - ^F 660 Plit 


O:f660 - O',' 


(A.9) 


(A. 10) 


With both parameters determined as a function of known val¬ 
ues, we can go back to Eq.|A.5|to obtain 


Hu+[NII] 




which is Eq.j^from Sect. 3.2 
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